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1. Introduction 

In random matrix theory, supersymmetry is an indispensable tool [H El El Hj . Recently, 
this method was extended from Gaussian probability densities to arbitrary rotation 
invariant ones. Presently, there are two approaches referred as superbosonization. 
The first approach is a generalization of the Hubbard-Stratonovich transformation for 
rotation invariant random matrix ensembles [5]. The basic idea is the introduction 
of a proper Dirac-distribution in superspace, extending earlier work in the context of 
scattering theory [UJ, universality considerations [7J, field theory [HI [Hj and quantum 
chromodynamics |10j . The second approach is the superbosonization formula developed 
in Refs. [TT| [T2] . It is an identity for integrals over superfunctions on rectangular 
supermatrices which are rotation invariant under an ordinary group. 

Here, we further extend the generalized Hubbard-Stratonovich transformation to 
the orthogonal and the unitary symplectic symmetry class in a unifying way. To this 
end, we use an analog of the Sekiguchi differential operator for ordinary matrix Bessel- 
functions. We also aim at a presentation which is mathematically more sound than the 
one in Ref. [5]. 

The article is organized as follows. The problem is posed in Sec. EJ We give an 
outline of the calculation in Sec. [31 In Sec. HI we present the generalized Hubbard- 
Stratonovich transformation. In Sec. G3, we carry out the calculation for arbitrary 
ensembles as far as possible. Then, we restrict the computation to the three classical 
symmetry classes. We, thereby, extend the supersymmetric Ingham-Siegel integral [5]. 
In Sec. |6j we give a more compact expression of the generating function in terms of 
supermatrix Bessel-functions. We show that the generating function is independent 
of the chosen representation for the characteristic function. The one-point and higher 
correlation functions are expressed as eigenvalue integrals in Sec. [71 In the appendices, 
we present details of the calculations. 

2. Posing the problem 

We consider a sub- vector space OJIat of the hermitian iV x iV-matrices Herm(2,iV). 
Herm (/3, N) is the set of real orthogonal ((3 = 1), hermitian ((3 = 2) and quaternionic 
self-adjoint ((3 = 4) matrices and (3 is the Dyson-index. We use the complex 2x2 
dimensional matrix representation for quaternionic numbers H. The results can easily 
be extended to other representations of the quaternionic field. For the relation between 
the single representations, we refer to a work by Jiang |13j . 

The object of interest is an arbitrary sufficiently integrable probability density P 
on OJItv- Later, we assume that P is an invariant function under the action of the group 




0(N) ,0 = 1 
U(AQ ,(3 = 2 
USp(2iV) , (3 = 4 



(2.1) 
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and 9J1 72 at = Herm (/?, N). Here, we introduce 72 = 1 for /3 e {1, 2} and 72 = 2 for /3 = 4 
and, furthermore, 71 = 2 / ~f 2 /l3 and 7 = 7172- These constants will play an important 
role in the sequel. 

We are interested in the fc-point correlation functions 

R k (x) = d k J P(H) Y[ tr 6(x p l N - H)d[H] (2.2) 



P =i 



with the k energies x = diag {x\, . . . ,x k ). Here, d is the inverse averaged eigenvalue 
degeneracy of an arbitrary matrix H G 9JIat. The measure d[H] is defined as in Ref. 
[14J , it is the product of all real and imaginary parts of the matrix entries. For example, 
we have d = 1/2 for O^tv = Herm (4, N) and d = 1 for no eigenvalue degeneracy as 
for m N = Herm(/?,iV) with (3 e {1,2}. We use in Eq. (Q the ^-distribution which 
is defined by the matrix Green's function. The definition of the /c-point correlation 
function (12.21) differs from Mehta's [15] . The two definitions can always be mapped onto 
each other as explained for example in Ref. [3]. 

We recall that it is convenient to consider the more general function 

R k (x (L) ) =d k I P(H) Y[tr[(x p + L p ie)t N - H^dlH] (2.3) 
J- p=i 



where we have suppressed the normalization constant. The quantities Lj in x^ = 
diag (xi + L%i£, ...,Xk + L k ie) are elements in {±1}. We define x ± = diag (xi ± 
is, . . . , Xk ± ie). Considering the Fourier transformation of ( 12.2ft we have 

r k (t) = (2n)- k ' 2 I R k (x)Y[exp(tx p t p )d[x] = 

Zu P=l 



p 

k « k 



= J P(H)f[txexp(iHt p )d[H} . (2.4) 

The Fourier transformation of (12.31) yields 

r k (t) = (2vr)- fc / 2 / R k (x^) J] exp (ix p t p ) d[x] = 

k 

= Yl l-L P 2mQ{-L p t p ) exp (eL p t p )} r k {t) (2.5) 
P =i 

where B is the Heavyside-distribution. 

As in Ref. [5], the fc-point correlation function is completely determined by Eq. 
( 12. 3p with L p = — 1 for all p if the Fourier transform (12.41) is entire in all entries, i.e. 
analytic in all entries with infinite radius of convergence. We obtain such a Fourier 
transform if the /c-point correlation function R k is a Schwartz-function on M. k with the 
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property 

/ \Rk{x)\ Y[exp (5x p ^ d[x] <oo , \/5 G R . (2.6) 

This set of functions is dense in the set of Schwartz-functions on R fc without this 
property. The notion dense refers to uniform convergence. This is true since every 

Schwartz-function times a Gaussian distribution exp I — e J , e > 0, is a Schwartz- 

es p=i / 

function and fulfils Eq. (12.61) . We proof that r k , see Eq. (12 .4p . is indeed entire in all 
entries for such /c-point correlation functions. To this end, we consider the function 

Tks{t) = J Rk{x) JJexp (iXptp) d[x], (2.7) 



p=l 



where 03^ is the closed fc-dimensional real ball with radius 5 G IR + . Due to the Paley- 
Wiener theorem [16], r k s is for all 5 G M + entire analytic. Let be another k- 
dimensional complex ball with radius S G M + . Then, we have 

lim sup \r k5 (t) - r k (t)\ < lim / \R k (x)\ TTexp ( 8x p ) d[x] = . (2.8) 

5->oo te rgc 8-*oo J \ J 

s 



The limit of r k s to r k is uniform on every compact support on C k . Thus, r k is entire 
analytic. 

The modified correlation function R k for all choices of the L p can be reconstructed 
by Eq. (I2.5p . In Sec. El we extend the results by a limit-value-process in a local convex 
way to non-analytic functions. 

We derive R k (x~) from the generating function 

f A det\H - (x~ + J V )1 N ] 

Zk + j) = j P (h) n lw fo n (") 

by differentiation with respect to the source variables [T7] 



n^sjp ,-> ' (210) 

where x _ + J = x _ <g> 1 4 + diag(Ji, . . . , J k ) ® diag (— 12, 12)- By definition, Z k is 
normalized to unity at J = 0. 



3. Sketch of our approach 

To provide a guideline through the detailed presentation to follow in the ensuing 
Sections, we briefly sketch the main ideas as in Ref. [3] and as further extended in 
the present contribution. 

To express the generating function (12. 9p as an integral in superspace, we write 
the determinants as Gaussian integrals over vectors of ordinary and Grassmann 



Random matrix ensembles and supersymmetry 



5 



variables. We then perform the ensemble average which is equivalent to calculating 
the characteristic function 

$(K) = J P(H)exp(itTHK)d[H} (3.1) 

of the probability density. The rotation invariance of P{H) carries over to §>(K). The 
ordinary matrix K contains the abovementioned vectors of ordinary and Grassmann 
variables as dyadic matrices. It has a dual matrix B in superspace whose entries are all 
scalarproducts of these vectors. The reduction in the degrees of freedom is fully encoded 
in this duality, as the dimensions of K and B scale with N and k, respectively. The 
crucial identity 

trK m = StiB m , Vm e N, (3.2) 

yields the super symmetric extension of the rotation invariant characteristic function, 

§(K) = $(tr K, tr K 2 , ...) = $(Str B, Str B 2 , ...) = , (3.3) 

which is now viewed as a function in ordinary and superspace. We rewrite it by inserting 
a proper Dirac-distribution in superspace, 

$(B)= J ®{p)5{p - B)d[p] (3.4) 

~ J J $(p) expfxStr (p - B)a]d[p]d[a] , (3.5) 

where the supermatrix p and a are introduced as integration variables. The vectors 
of ordinary and Grassmann variables now appear as in the conventional Hubbard- 
Stratonovich transformation and can hence be integrated out in the same way. We 
are left with the integrals over p and a. If we do the integral over p we arrive at the 
result 

Z k (x- + J) ~ J Q(a)Sdet ~ N ^(a - x~ - J)d[a]. (3.6) 

for the generating function. The superfunction Q is the superspace Fourier transform 
of $ and plays the role of a probability density in superspace, 

Q{a) = J $(p) exp(2Str pa)d[p] . (3.7) 

If we choose to integrate over a instead, we obtain another representation of the 
generating function 

Z k (x- + J) ~ J ®(p)I(p) exp[-zStr p(x~ + J)]d[p] , (3.8) 

which still contains the characteristic function. The distribution I(p) appears. It is the 
supersymmetric version of the Ingham-Siegel integral. It is a rotation invariant function 
resulting from the Fourier transformation of the superdeterminant in Eq. (13.61) . 

One way to proceed further is to diagonalize the supermatrix p and to integrate 
over the angles. We may omit Efetov-Wegner terms and have 

Z k {x- + J) ~ / ${r)I{r)(p{-vr,x- + J)d[r], (3.9) 
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where (p is a supermatrix Bessel-function. The differentiation with respect to J gives 
R k . We can introduce other signatures of L by Fourier transformation of Eq. (13.81) and 
identification with Eq. ( 12. 5ft . Eventually, we find the correlation functions R k - 



4. Generalized Hubbard— Stratonovich transformation 

In Sec. 14.11 we express the determinants in Eq. (12. 9p as Gaussian integrals and introduce 
the characteristic function of the matrix ensemble. In Sec. 14.21 we qualitatively present 
the duality between ordinary and superspace which is quantitatively discussed in Sec. 
14.31 Then, we restrict the matrix ensembles to the classical symmetry classes. In 
Sec. 14.41 we investigate the diagonalization of the dyadic matrix K appearing from the 
Gaussian integrals. The ambiguity of the supersymmetric extension of the characteristic 
function is discussed in Sec. 14.51 In Sec. 14.61 we present the symmetries of the appearing 
supermatrices. In Sec. 14.71 we replace the dyadic supermatrix in the supersymmetric 
extended characteristic function with a symmetric supermatrix discussed in the section 
before. 



4-1. Average over the ensemble and the characteristic function 

To formulate the generating function as a supersymmetric integral, we consider a 

2Nk 

complex Grassmann algebra A = ® Aj with iVA;-pairs {Qp, Qp}j,p °f Grassmann 

j=0 

variables [IB]. We define the k anticommuting vectors and their adjoint 

C P = (Ci„, • • • , Cn p ) T and (I = (Cr p , • • • , Cn p ) , (4-1) 

respectively. For integrations over Grassmann variables, we use the conventions of Ref. 
[T4"] . We also consider k iV-dimensional complex vectors {z p , ^}i< P <fe- In the usual 
way, we write the determinants as Gaussian integrals and find for Eq. ( 12.91) 

Z k (x- + J) = (-i) Nk J J d[(]d[z)d[H)P{H) x 

x exp ^ {Ctl H ~ ( x p + Jp^nKp + 4l H ~ ( x p ~ Jp)M^]^ ( 4 -2) 

k N k N 

where d[Q = U U d(j P d(* p , d[z] = U U dz jp dz* p and £ kN = C kN x A 2Nk . Using 
p=ij=i p=ij=i 

k 

{ClHCp + 4 Hz p) = tr H * (4-3) 

p =i 



with 



K = E M - CpCJ) (4.4) 



P =i 
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leads to 
Z k (x~ + J) 



TP (%(Wl N ;K)) x 



] d[Qd[z] . 



(4.5) 



p=i 



where the integration over H is the Fourier transformation of the probability density P, 



TP(n(Wl N ;K) 



P(H) exp UtvHK) d[H] 



(4.6) 



This Fourier transform is called characteristic function and is denoted by $ in Ref. [5] 
and in Eq. (13.11) . The projection operator 7t(9yt N ) onto the space £DT/v is crucial. For 
97t 72 7v = Herm iV) the projection operator is 

1 



^ Herm (J3,N);K 



2 I 



K + F(iT) 



with 



F(K) 

and the symplectic unit 

Y B = 



K 

(Y s <g> Ijv) K r (if <g> Ujv) 



/5 = 1 
= 2 
= 4 



(4.7) 



(4.i 



1 
-1 



(4.9) 



where Ijv is the N x iV-unit matrix. The transposition in Eq. ( 14. 8 j) can also be replaced 
by the complex conjugation due to 1ft = If. The projection onto the set of diagonal 



matrices @ K is 



7T 



AT 



diag JCn,^, . . . ,K 



NN 



(4.10) 



4-2. Duality between ordinary and superspace 

Is it always possible to find a supermatrix representation for the characteristic function 
TP such that Eq. (14.51) has an integral representation over supermatrices as it is known 
H2] for rotation invariant P on 9J1 72 tv = Herm (/?, iV)? The integral (14.51) is an 
integral over the supervectors Vj = (z* 1: . . . , z* k , —Qi, • • • , — (jk) T an d their adjoint 
Vj = (zji, . . . , Zjk, Cji, ■ ■ ■ , Cjk)- The adjoint "f" is the complex conjugation with the 
super symmetric transposition and "T" is the ordinary transposition. The entries of the 
matrix K are v^Vm. If we do not use any symmetry of the matrix ensemble, we can 
write these scalar products of supervectors as supertraces 

vly m = Strv m vl . (4.11) 
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Then, we can transform each of these supertraces with a Dirac-distribution to an 
integral over a (k + k) x {k + fc)-supermatrix. We defined the Dirac-distribution in 
superspace as in Refs. [T9| [ID]. The ambiguity discussed in Ref. [20] occurring by such 
a transformation is discussed in the subsections 14.51 and 16.31 

The procedure above is tedious. Using the symmetries of the ensemble (TP, Wl n), 
we can reduce the number of integrals in superspace. We will see that the number of 
commuting real integrals and of Grassmannian integrals is 2k 2 + 2k 2 (f3 = 2) or 4fc 2 + 4/c 2 
((3 G {1,4}) for a rotation invariant matrix ensembles on Herm(/3,iV). If there is not 
a symmetry the number of integrals has not been reduced. One has to integrate over 
N(N + 1) ordinary hermitian k x fc-matrices and their corresponding anticommuting 
parameters if the transformation above is used. 

4-3. Analysis of the duality between ordinary and superspace 

We consider an orthonormal basis {^Ik^ of WIn where d is the dimension of WIn- 
We use the trace tr A n A m = 5 nm as the scalar product and recall that 97tjv is a real 
vector space. Every element of this basis is represented as 

An = J2 X ^ne] n With Yl X % = 1 • ( 4 - 12 ) 

i=i j=i 
Here, e jn are the normalized eigenvectors of A n to the eigenvalues Xj n . Then, we 
construct every matrix H G 971 at in this basis 

d 

H = Y,KA n . (4.13) 

n=l 

We find for the characteristic function 



TP [% 



(m N ;K)^ = I p(^/i n A„]exp(\^MrAJf]4ff] 

J \n=l / V n=l / 



,n=l 



With help of Eq. (14.121) and an equation analogous to (14. lip , the characteristic function 
is 

TP (tt{WI n ; #)) = TP Str [J2 X inVe jn e) n V^ J A^j (4.15) 

with V = (f i, . . . , vn). We see that the matrix K is projected onto 

K = 7t(m N ;K) (4.16) 

where the projection is the argument of the characteristic function in Eq. (I4.14p . The 
matrices in the supertraces of (14.151) can be exchanged by (k + k) x (k + A;)-supermatrices 
with the Delta-distributions described above. If the ensemble has no symmetry then we 
have reduced the number of supermatrices to the dimension of 371^. Nevertheless, we can 
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find a more compact supersymmetric expression of the matrix K such that the number of 
the resulting integrals only depends on k but not on N. This is possible if K is a dyadic 
matrix of vectors where the number of vectors is independent of N and the probability 
distribution only depends on invariants of H. The ensembles with S0T 72 at = Herm (/3, N) 
and a probability density P invariant under the action of U ^' (N) fulfil these properties. 
It is known [51 [12] that these cases have a very compact supersymmetric expression. 
Furthermore, these ensembles are well analyzed for Gaussian-distributions with help of 
the Hubbard-Stratonovitch transformation [HEIE]. 

In the present context, the cases of interest are 9J1 12 n = Herm (/3, N) with a 
probability density P invariant under the action U^(iV). We need this symmetry to 
simplify Eq. ( 14. 151) . Let iV > 71 A;. This restriction also appears in the superbosonization 
formula [12]. If N < 71 k, one has to be modify the calculations below. For the 
superbosonization formula, Bunder, Efetov, Kravtsov, Yevtushenko, and Zirnbauer [20] 
presented such a modification. 

The symmetries of a function / carry over to its Fourier transform Tf. Thus, 
the characteristic function TP is invariant under the action of U ^' (N) . Let K Q 
be an arbitrary ordinary hermitian matrix in the Fourier transformation ( 14. 6 j) of the 
probability density. We assume that the characteristic function is analytic in the 
eigenvalues of K . Then, we expand TP as a power series in these eigenvalues. 
Since the characteristic function is rotation invariant every single polynomial in this 
power series of a homogeneous degree is permutation invariant. With help of the 
fundamental theorem of symmetric functions [21] we rewrite these polynomials in the 
basis of elementary polynomials. This is equivalent to writing these polynomials in 



the basis of the traces tr 



7T 



Herm (/?, N), Kq j , m E N. The analytic continuation 



of TP from Kq to K yields that the characteristic function in ( 14. 6 ft only depends on 



tr 



7r (Herm {{3,N),K 



Defining the matrix 







and its adjoint 



with 



V 



Y 



=•1) 



m e N. 



z k: Yz-y , . . . , Yz k , £1, . . . , Y^i . 



Yzi, . . . , Yzk, — Ci, 



1- 





Y T 



la/ 



/9 




1 

2 
4 



we find 



K = tt (Herm ((3,N);K 



7 



The crucial identity 

tr^vy 



str {yv^y 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



(4.21) 
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holds for all (3. It connects ordinary and superspace. For (3 — 2, a proof can be found 



in Ref. [SJ. In Appendix A we show that the equation 

StiViV 2 = StrV 2 V 1 
holds for all rectangular matrices of the form 



Vi 



At 



B 1 }c 
D 1 }d 



and Vo 



A 2 
C 2 



B 2 }a 
D 2 }b _ 



(4.22) 



(4.23) 



where Aj and Dj have commuting entries and Bj and Cj anticommuting ones. This 
implies in particular that Eq. (14.211) holds for all (3. Hence, we reduced the amount of 
supermatrices corresponding to K in Eq. (14.151) to one (2/c+2fc) x (2/c+2/c)-supermatrix. 
In Ref. [5], the characteristic function $ was, with help of Eq. (14.211) . extended 
to superspace. We follow this idea and, then, proceed with the Dirac-distribution 
mentioned above. 



4-4- Problems when diagonalizing K 

In Ref. [5] , two approaches of the duality relation between ordinary and superspace were 
presented. The first approach is the duality equation (14.211) for [3 = 2. In our article, 
we follow this idea. In the second approach, the matrix K was diagonalized. With the 
eigenvalues of K, a projection operator was constructed for the definition of a reduced 
probability density according to the probability density P. 

The latter approach fails because K is only diagonalizable if it has no degeneracy 
larger than j 2 . Moreover for diagonalizable K, one can not find an eigenvalue A = 0. 



This is included in the following statement which we derive in |Appendix E 



consisting of Grassmann variables r q = (jq , . . . , Tq ) T . Then, the matrix 



JV 

H {0) + J2l q [r q rl + Y(r^)] (4.24) 



Statement 4.1 

Let N,N e N, G Rerm((3,N), I e M.^ and {T q } 1<g< ^ ^ 2 N -dimensional vectors 

,\ 

// 

3=1 

can not be diagonalized H = f/diag (Ai, . . . , Xn)U* by a matrix U with the properties 

U^U = UU ] = 1 N , U* = Y{U) (4.25) 

and the body of U lies in U^(iV) iff has degeneracy larger than 72. Moreover, H 
has no eigenvalue A G M. 

In our particular case, K can not be diagonalized for k < N — 1. Hence, we do 
not follow the second approach of Ref. [5] . We emphasize that none of the other results 
in Ref. [3] is affected as they are proven by the correct first approach which we pursue 
here. 
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4-5. Ambiguity of the characteristic function in the supersymmetric extension 

In this section, we discuss the problem that the extension of the characteristic function 
TP from ordinary matrices to supermatrices is not unique. This results from the 
fact that symmetric supermatrices comprise two kinds of eigenvalues, i.e. bosonic and 
fermionic eigenvalues. Whereas ordinary symmetric matrices have only one kind of 
eigenvalues. In the supertraces, these two different kinds are differently weighted by a 
minus sign. To illustrate this problem, we also give a simple example. 

The rotation invariance of TP enables us to choose a representation TP® of TP 
acting on an arbitrary number of matrix invariants 

TP (tr K m \m e N) = TP(K) . (4.26) 

For this representation, a unique superfunction exists defined by 

$ ((7) = TP (Stra m |m e N) (4.27) 

where 

TP (Str £ m |m G N) = TP (tr K m \m G N) (4.28) 

with B = 7 _1 W t . However, the choice of the representation TP is not unique. The 
question arises whether it is a well defined object. It is clear that two representations 
TPq and TP\ are equal on Herm ((3, N) due to the Cayley-Hamilton theorem, 

TP (H) = TP^H) , H e Herm ((3, N). (4.29) 

The Cayley-Hamilton theorem states that there is a polynomial which is zero for H. 
Thus, H M with M > N is a polynomial in {H n }i< n < N . Plugging an arbitrary symmetric 
supermatrix a into the corresponding superfunctions $0 and $1 we realize that the 
choices are not independent such that 

$o(<t) t^i(<t) (4.30) 
holds for some a. 

For example with N — 2, k — 1 and (3 = 2, let the characteristic function 
TP(H) = TP (trH 3 ). We get with help of the Cayley-Hamilton theorem 

TPi (trH 2 , trH) = TP (2 tr H tr H 2 - tr 3 H) = TP (trH 3 ) = TP(H) . (4.31) 

Let the set of U < - /3 ^(p/g)-symmetric supermatrices be 

j<7 G Mat(7p/7g) a ] = a, a* = Y s (a) j and (4.32) 



Ys(o-) = { 



l2p 







Y s ® U, 



1 



(7 



(7 



l 2p 









Y} 



2q 



1 



1, 
2, 
4, 



(4.33) 



2q 



Random matrix ensembles and supersymmetry 



12 



with respect to the supergroups 

f U0Sp (+) (p/2g) , (3 = 1 
\]W(p/ q ) = i \J(p/q) , (3 = 2 . (4.34) 

[ U0Sp H (2p/ ? ) , (3 = 4 

Mat(7p/7g) is the set of {jp+jq) x (7p+7g)-supermatrices with the complex Grassmann 

8k 2 

algebra ® Aj. The definition of the two representations UOSp ^ of the supergroup 

UOSp can be found in Refs. [22j [H]. We refer to the classification of Riemannian 
symmetric superspaces by Zirnbauer [23J. 

We consider a U (l/l)-symmetric supermatrix a. This yields for the 
super symmetric extension of Eq. (14.311) 

TP Q (2Str aStr a 2 - Str V) + FP (Str a 3 ) = TP Q Q f 3 ^^" + Str 3(J ) ) ■ ( 4 - 35 ) 

One obtains the last equation with a theorem similar to the Cayley-Hamilton theorem. 
More specificly, there exists a unique polynomial equation of order two 

Str a 2 1 /„ 9 StrV 



o 2 --^— a- - StrV-^-f- ) =0 , (4.36) 
Str a 4 V StrV J V ' 

for a U (l/l)-symmetric supermatrix cr. 

The resulting integral in Sec. for the generating function Z k \<xn N =Herm(/3,N) is 
invariant under the choice of $o- This is proven in Sec. 16.31 Such an ambiguity of 
the supersymmetric extension of the characteristic function was also investigated by 
the authors of Ref. [20]. They avoided the question of the definition of a Dirac- 
distribution on superspace by the superbosonization formula. We introduce for the 
supersymmetric extension from Eq. (14.281) to Eq. (I4.27P a Dirac-distribution depending 
on the representation of the superfunction. 

4-6. Symmetries of the supermatrices 
We find for a chosen representation (FP$ 

Z k {x- + J) = {-i) k2N J $ (B) exp [-iStr (aT + J)B] d[(]d[z] . (4.37) 



(1*2 



N 



Here, we introduce ki = 72 A;, k\ = 71& and k = ^k. We will simplify the integral (14.371) 
to integrals over k\ eigenvalues in the Boson-Boson block and over fc 2 eigenvalues in the 
Fermion-Fermion block. 
For every (3, we have 

fi f = B , (4.38) 

i.e. B is self-adjoint. The complex conjugation yields 

YBY? , (3 e {1,4} 
YB*Y T , (3 = 2 
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with the (2k + 2k) x (2k + 2/c)-supermatrices 



Y 



(3=1 





l k 





lfc 











Y 



(3=4 



Y s <g> 1 A 
1* 
U fc 



(4.40) 



and Y 



(3=2 



diag (1, 0, 1, 0) (g> 1*. We notice that for the unitary case I? is effectively 

a (k + fc) x (A; + /c)-supermatrix, i.e. half the dimension. With help of the properties 
(14.381) and (I4.39[) we construct the supermatrix sets 



Eo(/3,fc) = < o e Mat(2fc/2fc) 



(7 



0", cr 




(4.41) 



A matrix in £ (/3, ^) fulfils the odd symmetry (I4.39p . We transform this symmetry with 
the unitary transformations 

" V2 l 2k 

U\(3=4 



u 



(3=1 



1 

7! 



1^ 

-it k 
o 



1^ 
«ifc 








V2 l 2 fc 



1 

7! 



o 







lfc 

-«lfc 





lfc 



(4.42) 



U\/3=2 = I4/C, according to the Dyson-index, arriving at the well-known symmetries 
of symmetric supermatrices [23], see also Eq. (14.321) . Defining the sets E (/?, k) = 
UTjq(/3, k)U\ we remark that the body of the Boson-Boson block of any element in 
these sets is a matrix in Herm ((3, k\). The body of the Fermion-Fermion block of any 
matrix in So(/3, k) lies in Herm (4//3, k 2 ). 

We introduce a generalized Wick-rotation e 1 ^ to guarantee the convergence of the 
supermatrix integrals. The usual choice of a Wick-rotation is e 1 ^ = 1 for investigations 
of Gaussian probability densities [3 HJ [2]. Here, general Wick-rotations [H] are also 
of interest. Probability densities which lead to superfunction as exp (— Stra 4 ) do not 
converge with the choice 1. Thus, we consider the modified sets 

Xl,(P,k) = %?: (P,k)%. (4.43) 

with ^ = diag (l 2 k, e"^ 2 ^)- Let k) be the set of supermatrices which contains 

only zero and first order terms in the Grassmann variables. 

In the sequel, we restrict our calculations to superfunctions which possess a Wick- 
rotation such that the integrals below are convergent. We have not further explored the 
set of superfunctions with this property, but we know that this set has to be very large 
and sufficient for our purposes. For example, superfunctions of the form 

$ ( a ) = 5(a) exp (-Str a 2n ) , n G N, (4.44) 

fulfil this property if ln<3>(cx) does not increase as fast as Strcr 2n at infinity. 
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Following Refs. [H El [TUj, $q(B) can be written as a convolution in the space of 
supermatrices £2,(/3,&) with a Dirac-distribution. We have 

Z k (x- + J) = {-i) k * N J J %(p)6 (p - UBrf) d[p] 

£k 2 N E$G8,fc) 



X 



x exp [— zStr 
where the measure is defined as 



x 



J)B] d[C]d[z 



d[p] = d[pi}d[p 2 } Yl drj nm d?i* nm . 



(4.45) 
(4.46) 



l<m<k2 
Kn<fci 



Here, {r] nm , r/* m } are pairs of generators of a Grassmann algebra, while p\ is the Boson- 
Boson and P2 is the Fermion-Fermion block without the phase of the Wick-rotation. 
Since p\ and p2 are in Herm (/?, ki) and Herm (4//3, fe), respectively, we use the real 
measures for d[pi] and d[p 2 ] which are defined in Ref. [H]. We exchange the Dirac- 
distribution by two Fourier transformations as in Refs. [5j [10]. Then, Eq. (14.451) 
becomes 



Z k (x~ + J) = (-t) k2N 2 2k ^ j j ^ (a)x 

x exp [zStr B (f/W - x~ - j)\ d[a]d{(\d[ 
where the Fourier transform of <J>o is 



z 



$o(p) ex P (— «Str pa) d[p] 



(4.47) 



(4.48) 



We write the supertrace in the exponent in Eq. (14.471) as a sum over expectation values 



1 N 

Str B (U ] oU - x- - J) = - tr *J ( uiaU ~ x ~ ~ J ) * 



7 r=i 



(4.49) 



with respect to the real, complex or quaternionic supervectors 

{Zjn, Zjni Cjn, Cjn} i< n < k 
\Zjn, 0, Cjn, 0}i< n <fc 



J 



°jn 



Z j+N,n 



— Z 



'j+N,n 

Z* 
/6 jn 



Cjn 
Cj+N,n 



~Cj+N,n 
/-* 



(3 



(4.50) 



Kn<k 



The integration over one of these supervectors yields 



I exp ~ tr (C/W — — J) ^ = ^ 2 Sdet " 1/7l p (<r - aT - J) 



(4.51) 



**2 
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p projects onto the non-zero matrix blocks of E_^,(/3, k) which are only (k + k) x (k + k)- 
supermatrices for [3 = 2. p is the identity for (3 G {1,4}. The Eq. (I4.5ip is true because 
U commutes with x~ + J. Then, Eq. (14.471) reads 



Z h {x- + J) = 2 2k ^ J ^$ (a)Sdet ~ N ^p (a - x~ - J) d[a] 



(4.52) 



£%(/3,fc) 



Indeed, this result coincides with Ref. [5] for (3 = 2 where the Fourier transform jF$ (o") 
was denoted by Q(o~). Eq. (14.521) reduces for Gaussian ensembles with arbitrary (3 to 
expressions as in Refs. [3] and [2j. The integral is well defined because e is greater 
than zero and the body of the eigenvalues of the Boson-Boson block is real. The 
representation (I4.52[) for the generating function can also be considered as a random 
matrix ensemble lying in the superspace. 

Eq. (14.521) is one reason why we called this integral transformation from the 
space over ordinary matrices to supermatrices as generalized Hubbard-Stratonovich 
transformation. If the probability density P is Gaussian then we can choose $o also 
as a Gaussian. Thus, this transformation above reduces to the ordinary Hubbard- 
Stratonovich transformation and the well-known result (I4.52p . 

5. The supersymmetric Ingham— Siegel integral 

We perform a Fourier transformation in superspace for the convolution integral (14.521) 
and find 

Z k (x- + J) = 2 2k ^ J ® {p)l!f> N \p)e W [-zStrp (x~ + j)] d[p] . (5.1) 

Here, we have to calculate the supersymmetric Ingham-Siegel integral 

I i k f3 ' N \p)= j exp(-iStipa + )Sdet- N/ ^pa + d[a] (5.2) 

with a + = a + zet^k- 

Ingham [2l] and Siegel [25] independently calculated a version of (15.21) for ordinary 
real symmetric matrices. The case of hermitian matrices was discussed in Ref. 
[26] . Since we were unable to find the ordinary Ingham-Siegel integral also for the 
quaternionic case, we give the result here. It is related to Selbergs integral [27]. Let 
R G Herm (/3,m), e > and a real number n > m — 1 + 2/(3, then we have 

J exp (-itiRS + ) det~ nhl S + d[S] = i- pmn/2 G { n P l mtm det x R S(R) (5.3) 

Herm (/3,m) 

where S + = S + zeIL 72m , the exponent is 
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and the constant is 



G 



09) 



n— m,m 



ry^ \ /3m(n—m+ 1)/2 
7T / 



n 



2tt^/ 2 

r (pj/2) 



(5.5) 



j=n— m+1 

r(.) is the Euler gamma-function and 0(.) is the Heavy side-function for matrices which 
is defined as 

1 , R is positive definite 
, else 



e(R) 



(5.6) 



The ordinary Ingham-Siegel integral was recently used in the context of supersymmetry 
by Fyodorov [26]. The integral was extended to the superspace E°/ 2 (2, k) in Ref. [5]. 
In this article, we need a generalization to all £%,(/3, fc), in particular f3 = 1,4. 

The integral (I5.2p is invariant under the action of U ^(ki/kz). Thus, it is convenient 



to consider I(r,e), where r = diag(r n , 



ri2, 



fc2; 



is the diagonal matrix of 



eigenvalues of p and contains nilpotent terms. The authors of Ref. [10] claimed in their 
proof of Theorem 1 in Chapter 6 that the diagonalization at this point of the calculation 
yields Efetov-Wegner terms. These terms do not appear in the p 2 integration because 
we do not change the integration variables, i.e. the integration measure d[p] remains the 
same. For the unitary case, see Ref. [5]. We consider the eigenvalues of p as functions 
of the Cartesian variables. We may certainly differentiate a function with respect to 
the eigenvalues if we keep track of how these differential operators are defined in the 
Cartesian representation. 

As worked out in Appendix C.l , the supersymmetric Ingham-Siegel integral (15.21) 
reads 

rW,N) 



/f ' W) (p) = CdetVie(ri)det fc r 3 exp (-e**etrr 2 ) 



l N 



6{r 2 ) 



|A fc2 (r 2 )|^ 



The constant is 
C 



-up 



7i 



k 2 N 



1_ 

'2tt 



M2 ('zl\ 



k-2 



\7i/ 



-) 



2k 2 {k 2 



(4//3) 
Vk 2 



with 



9k 2 



1 



r (2j/(3) 



^ ^(j-i)/^r (2/^) 

3=1 

while the exponent is given by 
N 72-71 

K = 

7i 2 
and the differential operator 

1 



(« 



A fe2 (r 2 ) 



det 



r N-b 
a2 



d 



dr, 



2 1 



u2 



h) ~^2 



(5.7) 



(5.J 



(5.9) 



(5.10) 



(5.11) 



l<a,6<fc 2 



is the analog to the Sekiguchi differential operator 



The complexity of D 



k 2 r 2 



b e) makes Eq. 



We derived it in Appendix B 
( 15.71) cumbersome, a better 



representation is desirable. To simplify Eq. (15.71) . we need the following statement 
which is shown in Appendix C.2 
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Statement 5.1 

We consider two functions F, f : Herm (4/(3, k 2 ) — > C invariant under the action of 

]j( 4 //3)^ 2 ^ anc ^ Schwartz-functions of the matrix eigenvalues. Let F and f have the 
relation 



F (p2) = f(p 2 )detp 2 



N/"fi—k 



for all p 2 e Herm (4//?, fc 2 ) 



(5.12) 



Then, we have 



F(r 2 )detV 2 |A fc2 (r 2 )| 4 / /3 exp (itrr 2 a 2 ) det^ 71 (e~^a 2 + ze^) d[a 2 ]d[r 2 ] 



k 2 Herm(4//3,fc 2 ) 



Wi/(0) 



F(r a )|Afa(r a )| 



4//3 



w 2 exp (ee 1 ^ tr r 2 ) -j^ 



<9 



AT-fei 



j2, 



d[r 2 ] (5.13) 



where the constants are 



/2vr 



w 1= - - 



VTi 



^ ^2fc 2 (fc 2 -l)//3 / N e -iil>N\ kl2 h. JL 



71/ 



W2 



^fclfca / 2vr \ 



ffc 2 

/., / \ 2fc 2 (fc 2 -l)//3 



nn 

6=1 a=l 



a 6—1 
— + 

7i 72 



yfc 2 



V7i/ \7i/ 



fc 2 fc 2 -i 

n 

J=0 



r(7v + i + 2j7/j) 
r(i + 2j//3) 



This statement yields for the supersymmetric Ingham-Siegel integral 

^2 / r> \ JV-fcl 



jf Ar) ( P ) = iye(r 1 )- 



where the constant reads 



jL 
2tt 



k\k2 



/2tt 
\7i 



A fc2 (r 2 )|^l = l \dr j2 

^ n 2fc 2 (fc 2 -l)//3 



<J(r 



J2J 



(5.14) 
(5.15) 

(5.16) 



- 



x^lln r(iv+i 



71/ 
- 2j/0) 



A" 



(JV-A;i)!7{ 



N 



fc 2 



-(4/0 

yfc 2 j=o 



r(i + 2j//?) 



(5.17) 



We further simplify this formula for (3 = 1 and (3 = 2. The powers of the Vandermonde- 
determinant A^(r 2 ) are polynomials of degree k 2 x 2(A; 2 — l)//3. The single power of 
one eigenvalue derivative must be 2(A; 2 — 1)/ (3 if we substitute these terms in Eq. (15.161) 
by partial derivatives of the eigenvalues, for details see | Appendix C.2j Hence, this power 
is a half-integer for (3 = 4. Also, A& 2 (r 2 ) has no symmetric term where all eigenvalues 
have the same power. Therefore, we can not simplify the quaternionic case in the same 
manner. 

We use the identities 



3=1 3 



n 9 , n A 4 (x) 

11 o 2(n-l) ' 
3=1 ° X j 



n-1 



n 



![(2n-2)rn(2j' + l) 

j=0 



(5.18) 
(5.19) 
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and find 



I^ N \p) = 2- k ^ 



_(N-2)\ 



x 



x 0(ri) detr\ 



(JV-l)/2 



n 



<9 



<9r 



iV-2 



and 



lf' N \p) = (_i) fe ( fe + 1 )/ 2 2- fc ( fc - 1 ) 



9(n) det rf J] 



,/2 



5{r 



J2> 



(5.20) 



x 



L(iV-l)!j 

AT-1 



9 



i2 



(5.21) 



For /3 = 4, we summarize the constants and have 



(N-k) 
x 0(ri) detrf +l ~ — 



4 fe A;! 



vr fc |A 2fc (r 2 )| 



2k 

n 



5 



J-2 



N-k 



S(r j2 ) . (5.22) 



These distributions are true for superfunctions whose Fermion-Fermion block 
dependence is as in Eq. (I5.12j) . Eqs. (I5.20p and (I5.2ip can be extended to distributions 
on arbitrary Schwartz-functions which is not the case for Eq. f)5.22p . The constants in 
Eqs. ( 15.201) and (15.211) must be the same due to the independence of the test -function. 

Statement 5.2 

Equations I15.2(J\) and h5.21\) are true for rotation invariant superfunctions $ which are 
Schwartz- functions in the Fermion-Fermion block entries along the Wick-rotated real 



axis. 



We derive this statement in Appendix C.3 



Indeed, the Eq. (I5.2ip is the same as the formula for the supersymmetric Ingham- 
Siegel integral for (3 = 2 in Ref. [5]. Comparing both results, the different definitions 
of the measures have to be taken into account. We also see the similarity to the 
superbosonization formula [HI El EE2J [HI [201 CD] for (3 G {1,2}. One can replace the 
partial derivative in Eq. (15.201) and (15.211) by contour integrals if the characteristic 
function $ is analytic. However for (3 = 4, more effort is needed. For our purposes, 
Eqs. (15.71) and (15.221) are sufficient for the quaternionic case. In the unitary case, the 
equivalence of Eq. (15.211) with the superbosonization formula was confirmed with help 
of Cauchy integrals by Basile and Akemann. [10J 



6. Final representation of the generating function and its independence of 
the choice for $ 

In Sec. 16.11 we present the generating function as a supersymmetric integral over 
eigenvalues and introduce the supersymmetric Bessel-functions. In Sec. 16.21 we revisit 
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the unitary case and point out certain properties of the generating function. Some 
of these properties, independence of the Wick-rotation and the choice of $o ; are a l so 
proven for the orthogonal and unitary-symplectic case in Sec. 16.31 



6.1. Eigenvalue integral representation 



The next step of the calculation of the generating function Zk(x + J) is the integration 
over the supergroup. The function & (p)ljf ,N \p) is invariant under the action of 

We define the supermatrix Bessel-function 



as in Refs. 



exp (Str sUrU^) dfi{U) 

U(«(fci/fc 2 ) 

. We choose the normalization 



(6.1) 



J f(a) exp (Str ax) d[e-^ /2 r ] }d[e^a 2 ]d[a 1 ] = 
J J /(a)^(s,x)|s^( ai ,e^ aa )|d[e^ 2 ]d[ Sl ]+b.t. 



(6.2) 



which holds for every rotation invariant function /. This normalization agrees with 
Refs. [BUJ EH ESI El Ej ■ The boundary terms (b.t.) referred to as Efetov-Wegner terms 
[321 1331 ITU] appear upon changing the integration variables [34J or, equivalently, upon 
partial integration [14J. The Berezinian is 

<(*i)<V^ 2 ) 



B 



08), 



sue^s 2 



V?(8ue^8 2 ) 



(6.3) 



fcl k2 

where Vk(si,e l ^s 2 ) = J] Yl ( s m — ^s m2 ) mixes bosonic and fermionic eigenvalues. 

n=l m=l 

These Berezinians have a determinantal structure 



Bf( Sl ,e^s 2 ) 



det 



det^ 



1 



Sal 



e^s h2 ) 2 



l<a<2fc 
Kb<k 



Sai - e^s b2 



J Ka,b<fc 



2 
4 



(6.4) 



For (3 = 2 this formula was derived in Ref. [32]. The other cases are derived in 
Appendix D We notice that this determinantal structure is similar to the determinantal 



structure of the ordinary Vandermonde-determinant raised to the powers 2 and 4. This 
structure was explicitly used [15j to calculate the fc-point correlation function of the 
GUE and the GSE. 
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We find for the generating function 
Z k (x~ + J) =2 2k{ - k -^ e ^ kl 1 1 rOW 




$ (r)/r* J,; (r) x 



x<p%> ka (-tr,x- + J) Bf{r u <*r 2 ) 



d[r 2 ]d[n] + b.t. . 



(6.5) 



The normalization of Z k is guaranteed by the Efetov-Wegner terms. When setting 
(k — I) with I < k of the source variables J p to zero then we have 

Z k (x- + J)\ Ji= =Jk=Q = Z^x- + 7) , (6.6) 

x = diag (xi, . . . jXi-i), J = diag (Ji, . . . , J/_i), by the integration theorems in Ref. 
[H ESI ESI EH El H3- This agrees with the definition flSfl . 

T/ie unitary case revisited 

To make contact with the discussion in Ref. [5], we revisit the unitary case using the 
insight developed here. 

For a further calculation we need the explicit structure of the supersymmetric 
matrix Bessel-functions. However, the knowledge of these functions is limited. Only 
for certain (3 and k we know the exact structure. In particular for (3 = 2 the 
supermatrix Bessel-function was first calculated in Ref. [32], [30] with help of the heat 
equation. Recently, this function was re-derived by integrating the Grassmann variables 
in Cartesian coordinates |14j . 

i k exp (— eStr r) 



<Pkk(- ir > x + J ) 



2 fc V 



det [exp (-ir mX {x n - J n ))] 1 



<m,n<k 



det [exp (ie^r m2 {x n + J n ))] x < 



<k 



(6.7) 



Si 2) (ri,e^r 2 )sj 2) (x — J,x + J) 
with x ± J = diag (xi ± J\, . . . , x k ± </&) and the positive square root of the Berezinian 

= s _ lV2 A fc ( Sl )A fc (e*s 2 ) 

r a i ~e^r b2 \ x< b<k 



Sf' 2) (n,e^r 2 ) = det 



_ 1 \fc(fe-l)/2 



(6- 



Vfc(si,e^s 2 ) 

(2) ('2') 

Due to the structure of <^, fc and -B^ , we write the generating function for p = 2 as an 
integral over $o times a determinant [5] 

1 



Z fc (x + J) 



where r r 



: _ l)fe (fe+l)/2 det - 



x b 



Ja — Jh 



J Ka,b<k 




$ (r) X 



x det[^j V (f mn ,x mn )©(r nil )exp(-£Strf TO „)] 1 < m)n < fc d[r 2 ]d[ri] + b.t. (6.9) 
diag (r 



['"mij ^ ^^2) 5 ^ti« diag (x m t/m; -^n "I - Jn) snd 

~ ,~ ^ , - - ~ y ° ^ 



8r m i exp ( iStr f mn x mn ) 
(N - l)\(r ml - e^r n2 ) 



—itp 







dr 



n2 



5{r n2 ) . (6.10) 
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Then, the modified fc-point correlation function is 



R k (x~) = / / $ (r) x 




x det [$ N (r mn , x mn )Q(r ml ) exp (-eStr r mn )} 
and the fc-point correlation function is 



l<m,n<k 



d[r 2 ]d[ri] +b.t(6.11) 




d[r 2 ]d[ri] +b.t. . (6.12) 



We defined x mn = diag (x m , x n ). The boundary terms comprise the lower correlation 
functions. The /c-point correlation function for (3 = 2 is a determinant of the 
fundamental function 



R R 

if there is one characteristic function J-'Pq with a supersymmetric extension $o 
factorizing for diagonal supermatrices, 



with $o : C — > C. For example, the shifted Gaussian ensemble in App. F of Ref. [5] is 
of such a type. 

In Eq. (16.131) we notice that this expression is independent of the generalized 
Wick-rotation. Every derivative of the fermionic eigenvalue r 2 contains the inverse 
Wick-rotation as a prefactor. Moreover, the Wick-rotation in the functions are only 
prefactors of r 2 . Thus, an integration over the fermionic eigenvalues r 2 in Eq. (16.111) 
cancels the Wick-rotation out by using the Dirac-distribution. Also, this integration 
shows that every representation of the characteristic function gives the same result, see 
Theorem 16.11 in the next subsection. However, the determinantal structure with the 
fundamental function in Eq. (16. 13j) depends on a special choice of $o- 

6.3. Independence statement 

For (3 = 1 and (3 = 4 we do not know the ordinary matrix Bessel-function explicitly. 
Hence, we can not give such a compact expression as in the case (3 = 2. On the other 
hand, we can derive the independence of the Wick-rotation and of the $o choice of the 
generating function. 

Statement 6.1 

The generating function Z k is independent of the Wick-rotation and of the choice of the 
characteristic functions supersymmetric extension <3>o corresponding to a certain matrix 
ensemble (P, Herm (/3, N)). 

Derivation: 

We split the derivation in two parts. The first part regards the Wick-rotation and the 
second part yields the independence of the choice of $o- 




(6.13) 



$ (r) = Sdet diag %{r u ), $ (r fcl ), $ (e # n 2 ) (e # r fc2 ) 



(6.14) 



5), 
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Due to the normalization of the supermatrix Bessel-function (16.21) . <f k ^ k (— ir, x~ + 
J) only depends on e"^. The same is true for <3>o- Due to the property 

(^lie) = e^D^% 2 (W) , (6.15) 
the Ingham-Siegel integral in the form (15. 7\i times the phase e l ^ kl ~ k2 ^ only depends on 
e"V 2 and e~ 1 ^ d / drz- The additional phase comes from the p-integration. Thus, we see 
the independence of the Wick-rotation because of the same reason as in the (3 = 2 case. 

Let $o an d $i be two different supersymmetric extensions of the characteristic 
function TP. Then these two superfunctions only depend on the invariants 
{Str cr mj }i<j<i and {Str c^li^i^u m j, n jJo,h £ N. We consider $ and $i as 
functions of C l ° — > C and C' 1 — > C, respectively. Defining the function 

A$(xi, . . .,x M ) = $ Q (x mi , . . . ,x miQ ) - $i(x ni , . . . ,x nii ), (6.16) 
where M = max{m a ,n(,}, we notice with the discussion in Sec. 14.51 that 

A$(a; 1 ,...,XM)U J =t r ^ = (6.17) 
for every hermitian matrix H. However, there could be a symmetric supermatrix o with 

A$(x 1 ,...,x M )\ Xj=StIa ^0. (6.18) 
With the differential operator 

k ^ {ie ll£) \ V k {r u e*n) ' { ] 

we consider the difference of the generating functions 
AZ k (x~ + J) = Z k (x~ + J)|$ - Z k (x~ + J)| $1 = 

= / lAfeCrOl^detVieCrOSrA^C^U^Ste^l d[n] (6.20) 

Here, we omit the Efetov-Wegner terms. The differential operator is invariant under the 
action of the permutation group Sftz) on the fermionic block Herm (4//3, k-z). Hence, 
we find 

M fry 

S) r A$(x)U. = strr4 2 =0 = d a(r)]J-Q-a-A$(x)\ Xj=Stlr3 

ae{0,...,N-ki} M i=l ■? 

|a|<fe 2 (iV-fei) 

M ffh 

d a ( ri )Y[-^-a-A<S>(x)\ Xj=tir3 = 

ae{0,...,V-fci} M j=l 3 

\a\<k 2 (N-ki) 

= 0, (6.21) 

where d a are certain symmetric functions depending on the eigenvalues r. At r 2 = 
these functions are well-defined since the supermatrix Bessel-functions and the term 
V A T 1 (r 1 , e"V 2 ) are C°° at this point. Thus, we find that 

AZ k (x~ + J) = 0. (6.22) 

This means that the generating function is independent of the supersymmetric extension 
of the characteristic function. □ 



T-2=0 
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7. One— point and higher order correlation functions 

We need an explicit expression or some properties of the supermatrix Bessel-function 
to simplify the integral for the generating function. For k = 1 we know the supermatrix 
Bessel-functions for all (3. The simplest case is (3 = 2 where we take the formula (16.121) 
with k = 1 and obtain 

Rl ( x ) = R^)( x , x ) = J J $ ( r ) M^ 2) dr 2 d n . (7.1) 

R R 

Since the Efetov-Wegner term in the generating function is just unity there are no 
boundary terms in the level density. For (3 £ {1,4} we use the supermatrix Bessel- 
function [221 ESI Q3| 



^21 ( — ir i x + J) = ex P I" - tStrr(x + J 

71 



2J 



x 



x [iStvr + J{r n -e^r 2 ) (r 21 - e^r 2 )] . (7.2) 

We find 

R 1 ( x -)=-i[ /$ (r)detri iV " 1)/2 Strr l r n- r 2il 



m — e*^r 2 ) 2 (r 2 i — e l ^r 2 ) 2 



1 / d \ N ~ 2 

x exp (-ixStrr) 6(n) — n I -^^r J S(r 2 )d[ ri ]dr 2 (7.3) 



for (3 = 1 and 

(ri — e^ri2) 2 (ri — e^r 22 ) 2 




R l{x -) = -4, / / * o(r)r r +1 st r r— -4 s ^_f^_ x 



x exp (-.-St,,-) e(r,)^ (4e-<)" «^<WH) 
for = 4. The differential operator has the explicit form 

* * V ( , 5) 



2 ' r2 dr 12 dr 22 2 r u - r 22 \dr 12 dr 22 / 
For the level density we have 

i?i(x) = / / $ (r) detr( iV_1 ^ 2 exp (— zsStrr) Strr- [ n o/ 21 ■ ; — — x 

KJ 2itJJ 0KJ 1 PV ; (r n -e^r 2 ) 2 (r 21 -e«/V 2 ) 2 

R 2 R 

1 / a \ ^-2 

1 / _../, a 



x (9(n) + 0(-n)) -^-^yy ^-^^rj 5(r 2 )d[n]dr 2 (7.6) 
for (3 = 1 and 

i? x (x) = -- / / $ (r)rr +1 exp(-zxStrr)Strr 7 ^'^^ Vt * 

7T y y (n - e^ri 2 ) 2 (ri - e^r 22 ) 2 

R R 2 

X (2iV + l)!r e U *») e ^r 12 -e^r 22 dNdri {7J) 
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for (3 = 4. The equations (17.4ft to (17.7ft comprise all level-densities for arbitrary matrix 
ensembles invariant under orthogonal and unitary-symplectic rotations. As probability 
densities which do not factorize are included, these results considerably extend those 
obtained by orthogonal polynomials. 

For higher order correlation functions we use the definition ( 12.3ft and the definition 
of the matrix Green's function. With help of the quantities L = diag (Li, . . . , G 
{±l} fc and L = L® l 2 j, this yields 

r r k ll P ' N) (lA exp f-eStr Lr 

R k (x)=2^ / / $ (r)hm E II L 



'e\0 ^— ' J-- 1 - J (27r«e _t ^ : 0* i 



X 



n-^iA(-^ (o) ^) 



(n,e^r 2 ) d[r 2 )d[ ri ) + b.t. (7.8) 

j=o 

for analytic correlation functions. We extend this formula to all rotation invariant 
ensembles by the universality of the integral kernel. First, we make a limit of a uniformly 
convergent series of Schwartz-functions analytic in the real components of its entries 
to every arbitrary Schwartz-function describing a matrix ensemble. The Schwartz- 
functions are dense in a weak sense in the sets of Lebesgue-integrable Functions L p and 
the tempered distributions. Thus, we integrate Eq. ( 17.8ft with an arbitrary Schwartz- 
function on IR fc and take the limit of a series of Schwartz-functions describing the 
ensembles to a tempered distribution which completes the extension. 



8. Remarks and conclusions 



We extended the method of the generalized Hubbard-Stratonovich transformation to 
arbitrary orthogonally and unitary-symplectically invariant random matrix ensembles. 
Due to a duality between ordinary and supersymmetric matrix spaces, the integral for 
the fc-point correlation function is over a superspace. This integral was reduced to an 
eigenvalue integral for all probability densities, including those which do not factorize. 
The results are in terms of the characteristic function. Thus, the characteristic function 
has to be calculated for the ensemble in question. Since the matrix Bessel-functions 
of the ordinary orthogonal and unitary-symplectic group [391 EH EQ] and, thus, the 
supermatrix Bessel-functions of UOSp (2k /2k) are not known explicitly beyond k = 1, 
we can not further simplify our results. However, we found the previously unknown 
determinantal structure of the Berezinian of UOSp (2k /2k). 

Up to the restriction N > ki, formula (17.81) is exact for every k, N and rotation 
invariant ensemble. Thus, it can serve not only as starting point for universality 
considerations [7] , but for all other studies. 

The expressions for the supersymmetric Ingham-Siegel integrals ( 15.201) . ( 15.211) and 
( I5.22p confirm the equivalence of the superbosonization formula [2TJI [TT] [T2"] with our 
derivation. A work for a proof of this equivalence for all /3's is in progress. The 
comparison of the superbosonization formula [T21 [TT] with Eq. (15.11) shows that the 
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crucial difference lies in the integration domain. However, the Dirac-distribution and 
the partial derivatives in the fermionic part imply a representation as a contour integral 
which is equivalent to the compact space used in the superbosonization formula. 
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Appendix A. Circularity of the supertrace for rectangular super mat rices 

The circularity for rectangular matrices of pure commuting entries or anticommuting 
entries was derived by Berezin [18]. Since we have not found the general theorem for 
arbitrary rectangular supermatrices, we give the trivial statement. 

Statement Appendix A.l 

Let the matrices V% and V 2 be the same as in Eq. \J^.23 ). Then, we have 



StiV 1 V 2 = StTV 2 V 1 (A.l) 

Derivation: 

We recall the circularity of the trace for rectangular matrices of commuting elements 
tr A\A 2 = tr A 2 A\ and its anticommuting analogue tr B\B 2 = — tvB 2 Bi which has been 
proven by Berezin [18]. We make the simple calculation 

StrViV 2 = trAaAj + tr-BiCs-trCi-Ba -tvDxD 2 
= tr A 2 A X - tr C 2 B 1 + tr B 2 d - tr D 2 D X 

= StrV 2 V! (A.2) 

□ 

For our purposes we must prove 

tr(UV) m = Str(VV j ) m . (A.3) 
We define V 1 = and V 2 = (VV^^V and get a = 2k, b = 2k, c = -? 2 N and 



d = 0. Applying corollary Appendix A.l| and reminding that tr A = Str A for a matrix 



of commuting elements and identification with the Boson-Boson block, we have the 
desired result flA.3[) . 
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Appendix B. A matrix— Bessel version of the Sekiguchi differential operator 

We derive a version for the Sekiguchi differential operator for the ordinary matrix Bessel- 
functions <p^p(y,x) on the connection between the Jack-polynomials and the ordinary 
matrix Bessel-functions. 

The Sekiguchi differential operator is defined as [28] 



D Ng (u,/3) = A^(z) det 
_d_ 



y N-b 
"a 



d 8 



Ajf(z) det 



- z, 



+ u) z, 



N-b 



2 

N-b 



l<a,b<N 

d 

+ u 

a 



J a r, 

az, 



(B.l) 



Ka,b<N 



Here, u is a boost and the expansion parameter to generate the elementary polynomials 
in the Cherednik operators, for more explicit information see Ref. |41j . Let J^\n,z) 
the Jack-polynomial with the partition n\ > ... > and the standard parameter 
a = i in Macdonald's [12] notation. The Jack-polynomials are eigenfunctions with 
respect to Dn z (u,/3) 

N 

D Nz (u,(3)J ( f(n,z)=H 



0=1 



n a + (N-a)^ + u 



(B.2) 



The aim is to find a similar differential operator for the ordinary matrix Bessel-function 
<PpP(y,x) such that 

N 



09) y_ 

N [ 12 



X 



,09) / ^ 



a=l 



x 



72 

09) f^. 
72' 



(B.3) 



Statement Appendix B.l 

The differential operator which fulfils Eq. W.3\) is 

d 



D% ) x (B) = A N 1 (x)det 



x 



.N-b 
a 



. (B.4) 



Ka,b<N 



Derivation: 

Kohler [43J has presented a connection between the Jack-polynomials and the matrix 
Bessel-functions. Let 



,2£ r , L 

z a = e l » and n a = —y a 

Z7T 



N + l 



— a 



2 



(B.5) 



then it is true 



Vn 



09) / f 



X 



lim 



72' / \A N (x)A N (y) 



A N (z) 



p/2 N 

Uz;K»-WjjP(n,z). (B.6) 

0=1 
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We expand the determinant in Eq. ( IB. II) and have 

D Nz (u,f3) = 



N 



,3 ( 8 
2 V'BZ + U 



e n 

me{0,l} N a=l 

Using the substitution (IB .51) and 
A(x) = 

we consider the limit 

f2m\ N 
lim — D Nz {u,f3) 

L— >oo \ Lt I 



N r 



o=l 



l—ma 



Zzsm [ — {x a — Xb) ) exp ( Z7T 

l<a<fe<iV 



(B.7) 



(B.8) 



1 



N 



lim 

L ^°°A(x) 



e n 



5/9 27TU 



2 \dx n 



A fx) x 



x 



N 

n 



A? 



l-m a 



e n 

P ( d 



P ( d 



2 \dx r 



+ tB 



1-2 
2 



A^(a;) det 



A" 1 ^) det 



x«-"i^- + (N-b)^-+tB 

ux a I X a / j l<a b<N 



_d_ 

8Xr 



iB 



l-m a 



Ka,b<N 



(B.9) 



Here, we defined a boost B = lim 2nu/L . The eigenvalue in Eq. ( IB. 21) is in the limit 



lim ( 



AT AT 



n 

o=l L 



n a + (N - a)~ + it 



(y a + 5) = det 1/72 <y + Bt 12N ) . (B.10) 



a=l 



We assume that Eq. (IB.6P is a uniformly convergent limit. Thus, we combine ( lB.6p . 
pi and (lB~T0l) with Eq. and find Eq. (jRij) . □ 

Indeed, D^l(B) is for the unitary case, (3 = 2, 



N f 8 \ 
D { l{B) = A N \x) J] f — + ,S J A^x) 

a=l ^ a ' 



(B.ll) 



Appendix C. Calculation of the supersymmetric Ingham— Siegel integral 



In Appendix C.I , we compute the Ingham-Siegel integral. We derive the statements 15. II 
and 15.21 in Appendix C.2 and Appendix C.3, respectively. 



Random matrix ensembles and supersymmetry 



28 



Appendix C.l. Decomposition of the Boson-Boson and Fermion-Fermion block 
integration 



We split a in its Boson-Fermion block structure 
pcr = 



e -#/2 



rr„ e""V 2 



(C.l) 



The following calculation must be understand in a weak sense. We first integrate 
over a conveniently integrable function and, then, perform the integral transformations. 
Hence, we understand ijf as a distribution where we must fix the underlying set 
of test-functions. For our purposes, we need Schwartz-functions analytic in the real 
independent variables. 

Since the superdeterminant of p (cr + lel^k) is 

det (<7i + letj.) 



Sdet pa" 



det 



e~^cr 2 + lett - 



a v (cri + ie±~ k ) 1 o\ 



(C.2) 



we shift cr 2 by analytic continuation to cr 2 + o n (o\ + islj:) 1 cr^ and obtain 



(p) 



y exp (-« tr ricri + 2 trr 2 o- 2 + 2 tr [r 2 a^ (cri + 



x 



£%(/3,fc) 

x exp (eStr r) 



det (e-^g-2 + teljfc) 
det (cri + lelfc) 



iV/71 



(C.3) 



An integration over the Grassmann variables yields 



(P) 



— 

"2V 



exp (eStrr) det r 2 x 



y exp (— itr riOi) det (cri + i£^' k ) 



d[(Ti] X 



Herm {fi,k\) 



J exp (ztr r 2 cr 2 ) det (e^cr 2 + lei*)* 771 ci[cr 2 ] 



(C.4) 



Herm (4//3,fc 2 ) 

With help of Eq. ( 15. 3ft we have 



r(/3,A0 



(P) 



-k 2 N n {P) 



Gjfi ( ) det K ri6(ri) exp (-e^etrra) 



2tt 



fcl&2 



x det fc r 2 y exp (2 tr r 2 cr 2 ) det Nhl (e"^cr 2 + zel s ) d[cr 2 ] 

Herm (4//3,fc 2 ) 

The remaining integral over the Fermion-Fermion block cr 2 , 

3T(r 2 ) = exp (-e^e tr r 2 ) y exp (i tr r 2 cr 2 ) det^ 71 (cr 2 + te^elj d[a 2 ] 



(C.5) 



(C.6) 



Herm (4//3,fc 2 ) 
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is up to a constant a differential operator with respect to r 2 times the Dirac-distribution 
of r 2 because the determinant term is for (3 G {1, 2} a polynomial in a 2 and for (3 = 4 
we use Cramers-degeneracy. We give several representations of this distribution. 

We first start with an eigenvalue-angle decomposition of o 2 = Us 2 U^ where s 2 is 
diagonal and U G U ^ A ^\k 2 ). Integrating over the group U < - 4///3 ^(fc 2 ), Eq. ( 1C.6I) becomes 

3{r 2 ) = exp (-e^etrr 2 )^ //3) x 

s 2 )det iV ^ (s 2 + ie^el~ k ) |A, 2 ( S2 )| 4 /^[ S2 ]. (C.7) 



For more information about the ordinary matrix Bessel-function 

¥>£ //J) (r 2 , s 2 ) = J exp (itrr 2 Us 2 tf) dfi(U) (C.8) 

U< 4 /«(fc 2 ) 

with normalized Haar-measure dfi(U) see in Ref. pSHig. The constant ^ is defined 
by 

J f(H)d[H] = g W I f(E)\A n (E)fd[E] (C.9) 

Hcrm (/3,n) M n 

independent of a sufficiently integrable function / which is invariant under the action of 
U ^\n). The Gaussian distribution is such a function. For the left hand side we obtain 

exp (-tr# 2 ) d[H] = 72 -«(a"-l)/2 2 -Mn-l)/4 7r n/2+/3n(n-l)/4 (Q 1Q) 
Herm (/3,n) 

The integral on the right hand side is equal to 



| exp (-7 2 |>n |A n (25)|M£]= < 

IB" > J =1 ' 



see Mehta's book [15]. Thus, we have 



r (Pj/2) 



2 -n(„-5)/4 j3 r (§ -h 1) 



2 -n(n-l)/2 7T n/2 T (j + 1) 
2 -n( 2 n-l/ 2 ) 7r n/ 2 T (2j + 1) , /? = 4 



= 1, 

/5 = 2, (C.ll) 



(C.12) 



This constant is the quotient of the volumes of the permutation group S(n) and of 
the flag manifold U^(n)/[U^(l)] n with the volume element defined as in Ref. [B] 
denoted by V0L3. 

We plug the differential operator of Appendix B (IB. 31) into Eq. (1C.7I) and have 
3(r 2 ) = g£ m exp (-e**etrr 2 ) {iji)~ k2N * 

r i N f 

<f («*7i*) / < //?) (r 2 , , 2 )|A fc2 ( S2 )| 4 ^[ S2 ] . (C.13) 



x 
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The integration over the eigenvalues leads to the Dirac-distribution 

*w=y U) V " 7l) 

and we find the representation for the supersymmetric Ingham-Siegel integral (15.71) . 
Appendix C.2. Derivation of statement s. 1\ 

The boost ie l ^e in the determinant can simply be shifted away because of 

D tIS («e # 7ie) exp (ee^trr 2 ) =exp (ee^trr 2 ) D%®{0) =exp (ee^trr a ) £>£{f (C.15) 

and Eq. (|C14B . Let & the set of U < - 4/ ' /3 - ) (A; 2 )-invariant Schwartz-functions on 
Herm (4/ /3, fca) — * C. The ordinary matrix Bessel-functions are complete and orthogonal 
in (5 with the sesquilinear scalar product 

(f\f) = [ f*(x)f'(x)\A k2 (x)\Wd[x) . (C.16) 



The completeness and the orthogonality are 



R-2 



k 2 

n si*, - x' p{j) ) 

fc fc 2 !^ 2) |A fe2 (x)|^|A fc2 (a;')|^ 1 ' ) 

where S(n) is the permutation group of n elements. We defined the constant 

S fi) W") ■ (ai8) 

Thus, we write D^Jf' in the Bessel-function basis 



x </ } y \<f>t m (y'WT\y')\ \^(y')\ A/ "d[y'] = 

= cf )_1 / det(i 7 iy) 1 ^^(y,x)^*(y, a 0|A fcl (y)| 4 //'d[ J/ ] (C.19) 
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with the action on a function / £ & 



x \A k2 (x')\ 4 ^\A k2 (y)\^d[x'}d[y} . (C.20) 

)(4/ 



Due to this representation of the Sekiguchi differential operator analog, i k2 D^^ is 



symmetric with respect to the scalar product (10. 161) 

{fll k 2D (m lf) = {t k 2D w) m {C21) 

Let L be a real number. Then, we easily see with help of Eq. ( IB. 41) 

<f det ^ = II ( L + | & - |) det x(L " lV71 • ( C - 22 ) 
Since the property ( 1C.21I) . we obtain for a function f E & 



det x L ^ \A k2 (x)\^D^f } f(x)d[x] 

R k 2 



= ("I)" 2 / /(x)|A fc2 (x)| 4 ^<fdetx^rf[ 3 

R fc 2 

= (-l)*jj^L + |&_:0 | /(^lA.^^l^detx^- 1 ^ 1 ^] . (C.23) 

6=1 R k 2 

The boundary terms of the partial integration do not appear because / is a Schwartz- 



function and D^/f^ has the representation ( 1C.19I) . 



Let F and / be the functions of statement 15.11 Then, we calculate 

J F(r 2 )det fc r 2 |A fe2 (r 2 )| 4 / /3 exp (itvr 2 a 2 ) det N ^ (e~^a 2 + iet~ k ) d[a 2 ]d[r 2 ] 

2 Herm (4//3,A;2) 

/(r 2 )det 7V/7l r 2 |A fe2 (r 2 )| 4 / /3 exp (*trr 2 a 2 ) x 



Rfe2 Herm(4//3,fc 2 ) 

xdet^ 71 (e-"V 2 + iel % ) d[a 2 ]d[r 2 ] 



Z ir) 9 ^ m J J /(r2)exp ^ trr2 ) i A ^)i 4//3 x 

R k 2 R k 2 

xdet^5 2 |A fc2 ( S2 )|^ (d£2> )' V 0if)(r2, 52 )rf[ S2 ]rfN = 

AT fc 2 



W« 



o=l 6=1 

X 



7i 72 



y y /(r 2 )exp (ee^trra) | A fe2 (r 2 )| 4 ^|A fc2 ( S2 )| 4 /^ 4 2 //3) (r 2 , S2 )rf[s 2 ]rf[r 2 ] 



K 2 M 2 



r d+2,/,3) /(0) ' (c - 24) 
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The second equality in Eq. (15. 131) is true because of 



k 2 



7 = 1 



d 



dr 



J2 



N-k-i 



f(r 2 ) exp (ee 1 ^ trr 2 ) det 



N/ji-k 



(C.25) 



r 2 =0 



The function in the bracket is F times the exponential term exp (ee 1 ^ tr r 2 J . 



Appendix C.3. Derivation of statement 1 5. S\ 
We have to show 

F(p 2 )det fc p 2 exp (iti p 2 a 2 ) &et Nhl a 2 d[a 2 }d[p 2 } 



Herm (4//3,fc 2 ) Herm (4//?,fc 2 ) 

k , ~ \N-2/0 

F ^U [-—) ^ ]2 )d[r 2 ] 

j=l \ 



J2 



(C.26) 



for every rotation invariant Schwartz-function F : Herm (4//3, fc 2 ) — > C and /3 G {1,2}. 
Due to 

y exp (ztrr 2 cr 2 ) deta^ 11 d[a 2 ] ~ ^ ^ y^exp [«rfc 22 tr(yl^ + v ^v)~\ d[v]dy x 

Herm(4//3,£! 2 ) K R4(fc 2 -l)/,3 

x y exp (i tr f 2 £ 2 ) detaf +2//3)/71 d[cr 2 ](C27) 

Herm (4//3,fc 2 -l) 

with the decompositions r 2 = diag (r 2 ,r fc22 ]l^) and 



a 2 



o 2 v 



we make a complete induction. Thus, we reduce the derivation to 




f(x)x kl y N exp [txtr(y + v^v)~\ d[v]dydx ~ j f(x) 

K R R 4(fc 2 -l)/,3 

where / : K — > C is a Schwartz-function. The function 

ki 



QN-2//3 

dx N ~ 2 /p 



f(y) = J f( x ) x 1 ex P ( tx v) dx 

is also a Schwartz-function. Hence, we compute 

f f(x)x hl y A r exp [«a;tr(2/ + t>^t>)] d[v}dydx 

R 4(fc 2 -l)//3 

/ [tr(j/ + v^v)] y N d[v]dy = 



(C.28) 



<J(x)d[x] (C.29) 



(C.30) 




(fc 2 -!)//3 



a v 2(fc 2 -l)//3 

[ -|J />(y + «M)d[i;]dy 



(fc 2 -l)//3 
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/(try)l^- 2 ^- 1 ^ 



AT-2(fc 2 -l)//3 



<5(x)<ix 



QN-2//3 



which is for /3 G {1,2} well-defined. 

Appendix D. Determinantal structure of the UOSp (2k/2k)— Berezinian 
Statement Appendix D.l 

Let k G N, Xi G C 2fc and x 2 G C fc . Xi and x 2 satisfy the condition 



(C.31) 



%al — %b2 



A 2fc (x!)At(x 2 ) 
V fc 2 (xi,x 2 ) 



(_l)fc(*-l)/2 det 



1 



X a l - Xf, 2 j l<a<2fc' [ (z al - X b2 ) 2 j l<<J<2fc 
K6<fe Kb<k. 



(D.l) 



(D.2) 



We prove this theorem by complete induction. 
Derivation: 

We rearrange the determinant by exchanging the columns 



det 



1 



Xal - %b2 J l<a<2fc 
Kb<k 



1 



(X a i -X 62 j Z J l<a<2k 
Kb<k . 



. l)fe (fe-l)/2 det 



_X a i — x o2 ' (x a i — x o2 )' 



l<a<2fc 



(D.3) 



Thus, the minus sign in Eq. (ID. 2p cancels out. 
We find for k = 1 

1 1 



det 



x n — x 2 (xn - x 2 ) 2 



(x n - X 2 i) 



. x 2 i - x 2 (x 2i - x 2 ) 2 . 



Xal - X fe2 (X a l-X 62 ) 



(xn - x 2 ) 2 (x 2 i - X 2 ) 2 



(D.4) 



. true. Let 






Si w 


\<a<2k 


V s 2 


Kb<k 





(D.5) 
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si 



s 2 



w 





r- 1 
1 
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%U ~ x 12 
1 





—^\2) 




. X21 — X12 


(x 2 l 


~ X 12 ) 2 
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1 




1 




_.T al — Xb2 








1 




1 




_X a l — Xu 

r 1 

1 


(Xal 


- X i2 ) 2 _ 

1 
1 




x ll — X b2 
1 


fall 


- ^62) 2 



3<a<2fc 
2<b<k 



and 



3<a<2fc 



2<b<k 



Then, we have 



det s = det Si det(s 2 — vs-t w 



X21 - Xb2 (x 2 i - x b2 ) 2 . 

(x u - x 2 i 



(DA) 



(x n - Xi 2 ) 2 (x 2 i - x 12 ) 2 
The matrix in the determinant is equal to 

(Xn - X a i) (X 21 - X a i) (x 12 - X b2 ) 2 1 



-1„..\T 



(s 2 - vs 1 w) 



(x al - Xi 2 ) 2 (xn - X o2 )(x 2 i - X b2 ) Xal ~ Xf, 2 

(X X1 - X al )(x 21 - X a i)(x 12 - X b2 ) P ab 
_ (X a i - Xi 2 ) 2 (xn - X 62 ) 2 (x 2 i - X 62 ) 2 (x a l ~ X b2 ) 



3<a<2k 
2<b<k 



where P ab is a polynomial 

Pab = [X a l — X b2 ){Xu — X b2 ){x\ 2 — X b2 ) — [X a i — Xi 2 )(xn — Xfe 2 )(x 2 i — x b2 ) 
-(X 2 l - X o2 )(x a l - X b2 )(xn - X i2 ) = 

= (Xn - x b2 )(x 2 i - x 62 )(xi 2 - X 62 ) + 

+ (x a l - X b2 ) [(Xn + X 2 l)(xi 2 + X 62 ) - 2xnX 2 i - 2xi 2 X o2 ] = 

= A i b l) + {x al -x b2 )Af ) . 



(D.6) 

(D.7) 
(D.8) 

(D.9) 



det(s 2 -vs^w) . (D.10) 



(D.11) 



(D.12) 



The polynomials A^ b and are independent of the index a. Due to the multilinearity 
and the skew symmetry of the determinant, the result is 

(xn - x 2 i) 



det s 



2k k 

II ( X H - X a l){.X 21 - Xal) EI ( X 12 - X b2 ) 
a=3 b=2 



(Xll-Xi 2 ) 2 (x 2 i-Xi 2 ) 2 2fc * 

11 {xai - x l2 y 11 (x n - x b2 y(x 21 - x b2 y 

a=3 6=2 



which completes the induction. 



dets 2 (D.13) 



□ 



Appendix E. Derivation of statement 14.11 

Let A be the wanted eigenvalue and is a commuting variable of the Grassmann algebra 
constructed from the {rq , Tq } Ptq . Then, we split this eigenvalue in its body A^ and 
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its soul i.e. A = X^ + A^. Let v the 7 2 iV-dimensional eigenvector of H such that 

Hv = Xv and v'v = 1 . (E-l) 

In this equation, we recognize in the lowest order of Grassmann variables that A^ -* is 
an eig envalue of # (0) . Then, let A (0) be an eigenvalue of the highest degeneracy 5 of 
lit ), i.e. S = dim ker(H^ — A^ljv). Without loss of generality, we assume that 
is diagonal and the eigenvalue X^ only appears in the upper left 5 x 5-matrix block, 

X^t s 



#(°> 

We also split the vectors in 6 and N — 5 dimensional vectors 







and r q = 






v 2 




T l2 



Thus, we find the two equations from ( IE. II) 
T xi v\ - A (1) wi + Ti 2 v 2 

T 2 lV! + 



v 2 








where T„ n 



N 

W 2 = 



_ X t N _ s + T 22 j 
+ ?fe T m )l- Eq. dESJ) yields 



(E.2) 

(E.3) 

(E.4) 
(E.5) 



'qn' qm 



H<® - Xl N _s + T> 



-22 



-1 



^21^1 . 



Hence, the body of t> 2 is zero and we have for Eq. ( IE. 41) 



Tuui - A (1) 



12 



# (0) - A1jv_ 5 + T, 



22 



-1 



T 21 ui = 



(E.6) 



(E.7) 



If the degeneracy is 5 > 72, we consider a 5-dimensional real vector w 7^ such 
that u^t>i = 0. Then, we get for the lowest order in the Grassmann variables of Eq. 
(IE.7I) times 



w 



(0) 



(E.8) 

where v x is the body of v\. The entries of vu^Tu are linearly independent. Thus, the 
body of v± is also zero. This violates the second property of ( IE. II) . 

Let the degeneracy 5 = 7 2 . Then, v\ is 7 2 -dimensional and is normalizable. For 
(3 = 4, we have the quaternionic case and the matrix before v± in Eq. (IE. 71) is a diagonal 
quaternion. Hence, it must be true 

A (1) l 72 = T u - T 12 [h® - Xt N _ s + T 22 ] ~* T 21 . (E.9) 

Considering the second order term in the Grassmann variables of Eq. ( IE. 91) . A's second 
order term is T n for (3 G {1,2} and trT n /2 for (3 = 4. Eq. (1E.9j) is unique solvable 
by recursive calculation. We plug the right hand side of Eq. ( 1E.9I) into the A^ on the 
same side and repeat this procedure. Hence, we define the operator 



0(ji) =— tr^Tn-r 12 
72 I 

O n+1 (f,) = 0[O n ^)] . 



H^-(X^+n)l N -s + T, 



22 



21 



and 



(E.10) 
(E.ll) 
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Then, A (1) = O n (A (1) ) is true for arbitrary n E N. The recursion is finished for n E N 
if A (1) = O n °(A (1) ) = O"°(0). Due to the Grassmann variables, this recursion procedure 
eventually terminates after the (7 2 iViV/2)'th time. Thus, the eigenvalue A depends on 
Grassmann variables and is not a real number. 
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